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Abstract 

As  a  model  for  streaming  multimedia  applications,  we  study  an  unreliable  retrial  queue  with 
infinite-capacity  orbit  and  normal  queue  for  which  the  retrial  rate  and  the  server  repair  rate  are 
controllable.  Customers  join  the  retrial  orbit  if  and  only  if  their  service  is  interrupted  by  a  server 
failure.  Interrupted  customers  do  not  rejoin  the  normal  queue  but  repeatedly  attempt  to  access 
the  server  at  i.i.d.  intervals  until  it  is  found  functioning  and  idle.  We  provide  stability  conditions, 
queue  length  distributions,  stochastic  decomposition  results,  and  performance  measures.  The 
joint  optimization  of  the  retrial  and  server  repair  rates  is  also  studied. 
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1  Introduction 


In  this  article,  we  analyze  an  M/G/1  retrial  queue  with  an  unreliable  server  whose  orbit  and 
normal  queue  have  infinite  storage  capacity  and  whose  retrial  and  server  repair  rates  are  control¬ 
lable.  Customers  in  service  join  the  retrial  orbit  if  and  only  if  they  are  interrupted  by  a  server 
breakdown  and  do  not  rejoin  the  normal  queue,  but  rather  attempt  to  access  the  server  directly  at 
random  intervals  independently  of  arrivals  or  other  retrial  customers.  However,  these  interrupted 
customers  can  regain  access  to  the  server  only  when  it  is  operational  and  idle  and  repeat  service 
until  they  have  been  completed.  Arriving  customers  who  find  a  failed  server  join  the  normal  queue. 
We  allow  for  both  active  breakdowns  which  occur  during  a  service  cycle,  and  idle  breakdowns  which 
occur  while  the  server  is  not  failed  but  idle.  The  server  may  not  breakdown  while  under  repair.  The 
times  between  customer  arrivals,  breakdowns,  retrials,  and  repairs  are  assumed  to  be  exponentially 
distributed  while  the  service  times  are  general. 

Over  the  past  two  decades,  advances  in  telecommunications  and  computer  networking  tech¬ 
nologies  have  reinvigorated  the  study  of  queueing  systems  in  general  and  retrial  queueing  systems 
in  particular.  The  model  we  present  here  is  well-suited  to  model  computer  network  streaming 
multimedia  applications.  The  primary  (or  normal)  queue  is  similar  to  a  1-persistent  carrier-sense 
multiple-access  (CSMA)  system.  When  the  oldest  packet  in  the  normal  queue  detects  that  the 
transmission  medium  (or  server)  is  free,  transmission  begins  immediately.  If  the  communication 
medium  fails  during  transmission,  the  packet  is  sent  to  a  retrial  queue  which  is  analogous  to  a 
non-persistent  CSMA  system.  If  the  medium  is  unavailable  (i.e.,  busy  or  failed),  then  the  retrial 
packet  waits  a  random  amount  of  time  before  checking  the  status  of  the  medium  again.  This  process 
repeats  until  the  retrial  packet  finds  the  transmission  medium  operational  and  idle.  In  this  sense, 
the  model  is  a  priority  queue  wherein  the  packets  that  are  not  interrupted  by  a  transmission  fail¬ 
ure  have  non-preemptive  priority  over  those  awaiting  availability  of  the  transmission  medium.  An 
important  application  is  that  of  streaming  voice  or  video  wherein  transmitted  packets  are  used  for 
playback  upon  reception  and  also  stored  for  future  use.  The  packets  used  for  immediate  playback 
are  time-sensitive  in  that,  if  they  are  not  received  within  a  given  time  threshold,  they  are  effectively 
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useless.  These  packets  correspond  to  the  priority  customers.  Packets  that  are  interrupted  during 
transmission  can  still  be  used  for  later  playback  from  the  stored  copy  of  the  stream,  but  their 
transmission  time  is  no  longer  important.  These  packets  constitute  customers  in  the  orbit.  By 
developing  analytical  expressions  for  congestion  and  delay  measures  in  stable  systems,  it  is  possible 
to  simultaneously  select  a  packet  retrial  rate  and  server  repair  rate  that  minimize  the  long-run 
average  cost  of  holding  customers  in  either  queue. 

In  addition  to  its  practical  relevance,  the  model  we  present  also  exhibits  extremely  interesting 
mathematical  properties  that  warrant  investigation  in  their  own  right.  The  presence  of  an  infinite 
waiting  space  for  primary  customers  introduces  an  interaction  between  the  two  infinite  queues. 
This  dynamic  does  not  exist  in  the  vast  majority  of  retrial  queueing  models  which  include  only 
an  infinite  retrial  orbit  and  do  not  consider  an  infinite  waiting  space  for  primary  arrivals.  We 
will  show  that  the  steady  state  orbit  size  and  the  overall  system  size  both  possess  a  stochastic 
decomposition  property.  Moreover,  an  interesting  stability  result  emerges,  namely  that  the  normal 
queue  may  remain  stable  even  if  the  condition  for  system  stability  is  violated.  Using  the  steady 
state  distributions  and  corresponding  queueing  performance  measures,  we  illustrate  the  means  by 
which  to  simultaneously  select  the  optimal  retrial  rate  and  repair  capacity  to  minimize  a  long-run 
average  cost  criterion. 

The  literature  addressing  retrial  queues  with  unreliable  servers  is  relatively  sparse  but  growing 
at  a  rapid  pace.  The  seminal  papers  in  this  area  are  [1]  and  [15].  All  models  considering  retrial 
queues  with  server  breakdowns  assume  an  M/G/1/1  loss  system  with  the  exception  of  [6]  and  [21]. 
Although  [10]  considered  an  M/G/1  retrial  queue  with  infinite-capacity  orbit  and  normal  queue, 
the  authors  did  not  consider  an  unreliable  server.  For  retrial  models  with  no  waiting  room  and 
server  breakdowns,  customers  arriving  to  find  the  server  unavailable  (busy  or  failed)  join  the  orbit. 
Some  models  (cf.  [2],  [3],  [7],  [16],  [20],  [23],  [25],  [27])  force  these  customers  into  the  orbit  while 
others  ([4],  [5],  [6],  [11],  [15],  [26])  provide  the  option  of  joining  the  orbit  or  departing  the  system. 
With  the  exception  of  two  cases  ([3]  and  [25]),  these  models  also  either  force,  or  provide  the  option 
for,  in-service  customers  interrupted  by  a  server  failure  to  join  the  orbit.  Our  model  differs  from 
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these  in  that  arriving  customers  who  find  the  server  busy  or  failed  join  the  normal  queue  whereas 
interrupted  customers  always  join  the  orbit  and  attempt  to  re-access  the  server  at  random  intervals. 
A  variety  of  failure  types  are  considered  in  the  literature  including  starting  failures  ([16],  [20],  [27]), 
vacations  ([7],  [23]),  active  breakdowns  ([6],  [25],  [26]),  and  like  our  model,  both  active  and  idle 
breakdowns  ([2],  [3],  [4],  [5],  [11],  [15]).  Most  orbits  are  assumed  to  behave  as  infinite-server  queues 
with  identical  exponential  service  times;  however  some  models  (cf.  [7],  [16],  and  [26])  consider 
orbits  as  FCFS  queues. 

For  retrial  systems  with  no  breakdowns  and  zero  capacity  in  the  normal  queue,  the  most 
common  optimal  control  strategies  include  the  optimal  routing  of  arriving  customers  ([9],  [10],  [22]) 
and  selection  of  the  optimal  retrial  rate  ([8],  [13],  [14]).  For  general  queueing  systems  (non-retrial 
queues),  researchers  such  as  [17],  [19],  and  [24]  have  considered  optimal  A^-policies  wherein  the  server 
remains  idle  until  exactly  N  (N  >  1)  customers  are  present  in  the  queue.  The  current  literature 
addressing  the  optimal  design  or  control  of  unreliable  retrial  queues  is  very  sparse.  It  appears  that 
only  [7]  formally  addressed  these  issues  for  a  retrial  queue  with  vacations.  In  that  work,  the  author 
presents  an  optimal  A^-policy,  an  optimal  T-policy  and  he  computes  the  optimal  retrial  rate  that 
minimizes  costs  using  an  A^-policy.  An  informal,  graphical  approach  to  the  optimal  control  and 
design  of  a  retrial  queue  with  vacations  was  presented  in  [20]  wherein  the  authors  examined  the 
impact  of  the  retrial  rate,  the  number  of  input  sources,  the  arrival  rate,  and  the  service  rate  on  the 
mean  waiting  time  and  throughput. 

As  a  model  for  streaming  multimedia  applications,  this  paper  is  concerned  with  the  analysis  and 
control  of  an  M/G/1  retrial  queue  with  an  infinite-capacity  orbit  and  normal  queue.  In  particular, 
we  consider  the  problem  of  simultaneously  selecting  an  optimal  retrial  rate  and  optimal  repair  rate 
with  the  objective  of  minimizing  the  long-run  average  operating  cost  which  is  a  function  of  the 
key  queueing  performance  measures.  To  this  end,  using  the  method  of  supplementary  variables 
and  a  classical  generating  function  approach,  we  derive  the  steady  state  joint  distribution  of  the 
orbit  size  and  normal  queue  size  when  the  server  is  idle  (operational  and  not  occupied),  failed 
(non-operational  and  being  repaired),  or  busy  (operational  and  occupied).  Using  these  results,  we 
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obtain  the  joint  generating  function  of  the  orbit  size  and  normal  queue  size  as  well  as  the  generating 
function  for  the  overall  system  size  (the  total  number  of  customers  in  orbit,  normal  queue  and  in 
service),  independent  of  the  server’s  status.  We  provide  a  necessary  and  sufficient  condition  for 
stability  of  the  orbit  and  system  as  well  as  a  (distinct)  condition  for  stability  of  the  normal  queue. 
Moreover,  we  show  that  the  steady  state  length  of  the  retrial  queue  and  the  system  size  can  be 
stochastically  decomposed. 

The  remainder  of  the  paper  is  organized  as  follows.  Section  2  provides  the  model  description 
and  mathematical  notation.  In  section  3  we  establish  stability  conditions  and,  by  means  of  gener¬ 
ating  functions,  derive  the  queue  length  distributions,  key  queueing  performance  measures,  as  well 
as  the  steady  state  distribution  of  the  server’s  status.  In  section  4  we  present  stochastic  decompos- 
ability  results  for  queue  length  distributions.  Finally,  section  5  presents  and  illustrates  a  nonlinear 
optimization  problem  for  the  optimal  selection  of  the  retrial  and  server  repair  rates. 

2  Model  Description 


Customers  arrive  to  the  system  according  to  a  homogeneous  Poisson  process  with  rate  A  >  0. 
Service  times  form  an  independent  and  identically  distributed  (i.i.d.)  sequence  of  random  variables 
with  absolutely  continuous  distribution  function  (d.f.)  B,  probability  density  function  (p.d.f.)  6, 
and  service  completion  rate 


/i(x) 


b{x) 

l-B{xy 


X  >  0. 


For  s  >  0,  let 


^  r 


%{x)dx 


Jo 

denote  the  Laplace  transform  of  b.  Server  failures  occur  according  to  a  Poisson  process  with  rate 
^  >  0  when  the  server  is  not  being  repaired.  The  repair  time  is  exponentially  distributed  with 
rate  parameter  a  >  0.  An  in-service  customer  interrupted  by  a  server  failure  enters  the  orbit  and 
spends  an  exponential  amount  of  time  there  with  rate  9  >  0,  after  which  it  either  enters  service 
(if  possible)  or  remains  in  the  orbit  for  an  additional  exponentially  distributed  time  with  rate  6. 
The  arrival,  service,  failure,  repair,  and  retrial  processes  are  assumed  to  be  mutually  independent. 
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Denote  by  Qt  the  number  of  customers  in  the  normal  queue  at  time  t,  excluding  any  customer  that 
might  be  in  service,  and  let  Rt  denote  the  number  of  customers  in  the  orbit  at  time  t.  The  random 
variable  Ut  is  the  occupation  status  of  the  server  given  by 


^  _  f  1,  if  the  server  is  occupied  at  time  t 

\  0,  if  the  server  is  not  occupied  at  time  t 

while  St  describes  the  operational  status  of  the  server  at  time  t  defined  by 

^  f  1,  if  the  server  is  not  failed  at  time  t 

*  1  0,  if  the  server  is  failed  at  time  t 


Let  Xt  denote  the  elapsed  service  time  of  the  customer  in  service  at  time  t  so  that  the  continuous¬ 
time  stochastic  process,  {{Qt,Ut,  Rt,  St,  Xt)  :  t  >  0}  describes  the  state  of  the  system.  Let  Nt 
denote  the  total  number  of  customers  in  the  system  at  time  t  (i.e.,  in  orbit,  normal  queue,  and  in 
service).  Assume  that  as  t  ^  oo,  Qt  ^  Q,  Rt  ^  R,  St  ^  S  and  Nt  ^  N  where  denotes  weak 
convergence. 

Now  define 


7ro,o,i,i  =  lim  P{Qt  =  0,Ut  =  0,  Rt  =  j,  St  =  1),  j  >0 

C— >oo 

vrfc,o,i,o  =  lim  P{Qt  =  k,Ut  =  0,  Rt  =  j,  St  =  0),  j,k>0 

C— >oo 

=  lim  P{Qt  =  k,Ut  =  I,  Rt  =  j,  St  =  I,  Xt  <  x),  j,k>0 

t^oo 

as  the  limiting  probabilities  that  the  system  is  in  an  idle,  failed,  or  busy  state,  respectively.  With 
the  transform  variables  zi  and  Z2  corresponding  to  the  orbit  size  and  normal  queue  size,  define 

OO 

00,0,1(2:1)  =  ^2:^710,0,7,1, 

7=0 

OO 

0fc, 0,0(21)  ^  ^  2]^7rfc,o,j,o, 

7=0 

OO 

4>k,l,l{x,Zl)  =  y~]2:j7rfc,l,j,l(x). 

7=0 

These  are,  respectively,  the  generating  functions  for  710,0,7,1,  7ifc,o,7,o,  and  TTk,i,j,i{x)  with  respect  to 
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the  orbit  size.  Further  define,  respectively, 

OO 

k=0 

OO 

^pl^l{x,Zl,Z2)  =  '^Z2(pk,l,l{x,Zl), 
k=0 

the  generating  functions  for  4'kfl,o{zi)  and  4>k,i,i{x,  zi)  with  respect  to  the  normal  queue  size.  The 
joint  p.g.f.  of  the  orbit  and  normal  queue  size  when  the  server  is  not  failed  and  busy,  is  given  by 

roo 

if^lAzi,  Z2)  =  /  'llJl^l{x,Zl,Z2)dx. 

Jo 

Let  p  denote  the  joint  probability  mass  function  (p.m.f.)  of  R  and  Q  while  q  denotes  the  p.m.f.  of 
N.  By  summing  over  the  three  distinct  and  exhaustive  server  states,  we  denote  by 

CXD  OO 

G{zi,Z2)  =  k)z{z2  =  (t)o,o,i{zi)  +  i^o,o{zi,Z2)  +  '4)l^l{zi,Z2), 

k=0  j=0 

the  joint  generating  function  for  the  orbit  and  normal  queue  size.  In  a  similar  manner,  we  denote 
by 

OO 

Hiz)  =  '^q{j)z^  =  (pofiAz)  +  'ipo,o{^,z)  +  ZTpi^i{z,  z), 
j=o 

the  generating  function  for  the  overall  system  size. 

In  the  next  section,  we  provide  stability  conditions  and  formally  derive  the  generating  functions 
defined  in  this  section.  Subsequently,  we  use  these  to  characterize  queue  length  distributions  and 
performance  measures. 

3  Stability  Analysis  and  Steady  State  Equations 


In  this  section,  we  provide  a  necessary  and  sufficient  condition  for  stability  of  the  overall  queue¬ 
ing  system  and  derive  the  steady  state  joint  distribution  of  the  orbit  and  normal  queue  size  when 
the  server  is  idle,  failed,  or  busy,  respectively.  Subsequently,  we  obtain  the  joint  distribution  of 
the  orbit  size  and  normal  queue  size,  and  the  distribution  of  the  system  size,  independent  of  the 
server’s  status.  Additionally,  we  obtain  standard  queueing  performance  measures  as  well  as  the 
limiting  distribution  of  the  server’s  status. 
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Before  proceeding  to  the  main  result,  we  first  provide  a  lemma  that  is  needed  to  characterize 
the  stability  conditions  and  steady  state  distributions.  As  in  Aissani  and  Artalejo  [4],  define  the 
fundamental  server  period  as  the  time  from  which  a  service  cycle  begins  until  the  next  time  at 
which  the  server  is  able  to  initiate  a  new  service  cycle.  Denote  this  random  duration  by  T.  Let  Nr 
and  Nq  respectively  denote  the  number  of  customers  entering  the  orbit  and  normal  queue  during 
(0,T],  and  let  a{i,j)  =  P{Nr  =  i,Nq  =  j),  i,j  >  0.  Define  the  generating  function 

OO  OO 

Q(zi,Z2}  =  l^il  <  1,  |2;2|  <  1. 

i=0  j=0 

Then  one  can  verify  (see  [4])  that 


Q{zi,Z2) 


b*{C  +  K^-Z2))  + 


aziCjl  -  b*{C  +  A(1  -  Z2)) 
{a  +  A(1  —  Z2)){^  +  A(1  —  Z2)) 


Now  let  us  define  the  quantity 


£  =  0 


A(l-y(0)(a  +  0 


where  h*{^)  is  the  Laplace  transform  of  the  service  time  p.d.f.  evaluated  at  Using  these  definitions 
we  have  the  following  important  result  which  is  needed  to  obtain  our  main  results. 


Lemma  1  For  either  \zi\  <  1  or  \zi\  <  1  and  pi  >  1  the  relation 

Z2  -  Q(zi,Z2) 

has,  as  a  funetion  of  Z2,  one  and  only  one  zero,  g{zi),  inside  the  region  \z2\  <  1.  In  ease  zi  =  1, 

g{l)  is  the  smallest  positive  real  zero  with  g{l)  <1  if  pi  >  1,  and  p(l)  =  1  i/pi  <  1. 

Proof.  The  proof  is  similar  to  that  of  Theorem  3  in  [18,  p.  351-352].  Now  for  the  first  part, 
on  \zi\  <  1,  applying  Rouche’s  theorem  to  the  function  Z2  and  the  generating  function  Q{zi,Z2), 
we  conclude  that  there  is  one  and  only  one  zero,  g{zi),  for  each  Z2  inside  the  unit  disk  |2;2|  <  1- 
For  the  second  part,  consider  the  quantity  pi  and  the  case  when  zi  =  1.  The  function  Q{zi,Z2)  is 
monotonically  increasing  in  Z2  for  Z2  G  [0, 1]  such  that  0  <  <5(1,0)  <  1  and  <5(1, 1)  =  1.  Thus,  if 
zi  =  1  and  Pi  >  1,  then  g(l)  is  the  minimal,  positive  real  zero  with  g(l)  <  1.  On  the  other  hand, 

if  pi  <  1,  5(1)  =  1  is  the  unique  zero,  and  this  completes  the  proof.  ■ 
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Using  Lemma  1,  we  now  characterize  the  stability  condition  for  the  overall  system  size  (and  orbit 
size),  as  well  as  the  generating  functions  V’o,o(-2i) -22),  and  2:2).  Theorem  1  states 

that  the  long-run  proportion  of  time  the  server  is  available  for  serving  customers  must  exceed  the 
long-run  proportion  of  time  the  server  is  busy  if  the  system  is  to  remain  stable. 


Theorem  1  The  queueing  system  is  stable  if  and  only  if  p  <  1  where 

Kl-b-ma  +  O 

»= — — ■  w 

In  sueh  a  ease,  the  generating  funetions  <^0,0,1(21);  V’o,o(-2i,  2^2)?  o,nd  Z2)  are,  respeetively, 

given  by 

" « exp (-1  f  1  _  P) 

(C)(a  +  0  [  dJz,  9{u)-u  J 

,  ^  ,  f  igjzi)  -  Zi)[a  +  A(1  -  g(2:i))][^  -t  A(1  -  Z2)][z2  -  B{z2)  -  zijl  -  B{z2))] 

\  {Z2- B{z2))[a  +  \(1- Z2)][i  +  \{1- Z2)]- ai{l- B{z2))zi 

Xzi{l- B{z2)){z2  -  zi){l  -  g{zi))[a  +  ^  + X{1- g{zi))]  \  I.(f)o,o,i{zi) 


{Z2  —  B{z2))  [a  +  A(1  —  Z2)]  +  A(1  —  Z2)]  —  0^5(1  —  B{z2))zi  j  (9(^1)  —  Zi)[a  +  A(1  —  ^(2:1))] 


(3) 


and 


1pl,l{zi,Z2)  = 


{Z2  -  Zi){l  -  g{zi))[a  -I-  C  -I-  A(1  -  g{zi))][a  +  A(1  -  Z2)] 


{z2  -  B{z2))  [a  +  A(1  -  Z2)]  [C  +  A(1  -  Z2)]  -  0^(1  -  B{z2))zi 

(1  -  Z2){g{zi)  -  zi)[a  +  A(1  -  g{zi))][a  C  +  A(1  -  Z2)]  \  A(1  -  B{z2))(l>o,o,i{zi) 

(22  —  B(z2))  [a  -I-  A(1  —  22)1  [C  +  A(1  —  22)]  —  0^(1  —  .B(22))2i  J  (5(21)  —  2i)[q;  -I-  A(1  —  3(21))] 


where 


B{z2)=b*{i  +  \{l-Z2)) 


and,  for  zi  G  [0, 1],  g{zi)  verifies 


g{zi)  =b*{^  +  A(1  -  g{zi)))  -h 


af,zi[l  -  b*{^  +  A(1  -  g{zi)))] 


[a  -h  A(1  -  g{zi))]  [(  +  A(1  -  g{zi))] 
Proof.  For  j  >  0  and  A:  >  1  with  Trk^i-ig{x)  =  0,  the  balance  equations  are 

roo 

{a  +  A)7ro,o,j,o  =  C72o,o,i,i  +  ^  /  7ro,ij-i,i(x)(ix, 

Jo 

roo 

(o  -|-  A)7r/j  Q  J  g  A7r^_]^  Qj  Q 

Jo 

roo 

(A  +  ^  +  J0)vro,oj,i  =  avrg^oj^o  -h  /  p{x)T:Qpjp{x)dx, 

Jo 

^TTgpj- i(x)  =  -[p{x)  -h  A  -F  C]7ro,lJ,l(x), 


(5) 

(6) 
(7) 


d^ 

dx 


72fc,lJ,l(x)  =  -[p{x)  -h  A  -h  C]7rfc,ij,i(x)  -h  XTTk-ipjpix) 


(8) 

(9) 


with  boundary  conditions 


7ro,ij,i(0)  =  avri,oj,o  +  A7ro,oj,i  +  /  /i(x)7ri4j,i(x)dx  +  (j  +  l)6'7ro,oj+i,i,  (10) 

Jo 


poo 

^fc,lj,l(0)  CK^fc+l,0j,0  • 

Jo 


(11) 


Multiplying  both  sides  of  Eqs.  (5)  through  (11)  by  and  summing  over  all  j,  we  obtain,  respec¬ 
tively,  the  following  equations: 


roo 

(a -h  A)(/)o,o,o(^i)  =  ^</>o,o,i(2:i) +^2^1  /  (l)o,i,i{x,zi)dx, 

Jo 

poo 

(a  +  A)(/>fc,o,o(2i)  =  A(/)fc_i, 0,0(21)  +  ?2i  /  (pk,i,i{x,  zi)dx, 

Jo 

d  f°° 

(A  -h  O'('o, 0,1(21)  -h  021— (/)o, 0,1  (21)  =  a(?!)o,o,o(2i)  -h  /  ^(x)(/)o,i,i(x,  2:i)(ix, 

«2i  Jo 


d 

— (/)o,i,i(x,  Zi)  =  -[//(x)  -h  A  -h  i](t)o,i,i{x,  Zi), 

JP 

dx 


(12) 

(13) 

(14) 

(15) 


4>k,i,i{x,  21)  =  -[/r(x)  -h  A  -h  ^](/>fc,i,i(x,  zi)  +  A0fc_i,i,i(x,  zi),  (16) 


with  boundary  conditions 

d 

0o,i,i(0, 21)  =  a(/)i, 0,0(21) -h  A(/)o, 0,1(21) -h  /  /u(x)(/)i,i,i(x,2:i)(ix-h0—(?!)o, 0,1(21),  (17) 

Jo  ClZi 

poo 

(i>k,i,i{0,zi)  =  a(j)k+i,o,o{zi)  +  n{x)(j)k+i,i,i{x,zi)dx.  (18) 

Jo 

Multiplying  both  sides  of  Eq.  (12)  by  (13)  by  and  summing  over  all  A:  >  0,  we  obtain 

poo 

[a  + X{1  -  Z2)]ipo,o{zi,Z2)  =  ((pofipzi)  +  Czi  ipiPx,zi,Z2)dx.  (19) 

Jo 

Performing  a  similar  operation  on  Eqs.  (15)  and  (16)  as  well  as  (17)  and  (18),  we  obtain,  respec¬ 
tively, 

fi 

(20) 


d 

—  •01,l(x,  Zl,  Z2)  =  -[px)  -h  ^  -h  A(1  -  Z2)]'lpl,l{x,  Zl,Z2) 


and 


d  ex 

l/:i,l(0,  Zi,  Z2)  =  A(/)o, 0,1(21)  -h  0^(/>O,O,l(2l)  H - [V'0,0(2l,  22)  -  (/>0,0,o(2l)] 

CLZl  Z2 


1 

H - /  ;u(x)[V’i,i(x,2:i,2:2)  -  <('o,i,i(a:,2i)](ix.  (21) 

22  Jo 
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Solving  Eq.  (20),  we  obtain 


Using  Eqs.  (14),  (19),  and  (22)  in  (21)  we  obtain 


'lpl^l{0,Zi,Z2)  = 


^(^2  -  -  [^(1  -  2:2)  +  ?(1  -  a+A(l-^2))]^0'04(2^l) 


2^2-  b*{^  + X{1  -  Z2))  + 


Q.^lg(l-b*(g+A(l-^2)) 

(q+A(1-2:2)){?+A(1-22)) 


Now  by  Lemma  1,  the  denominator  of  Eq.  (23)  has,  for  any  zi  in  the  unit  disk,  a  zero  in  the 
region  \z2\  <  1-  This  must  also  be  a  zero  for  the  numerator;  therefore,  we  have  from  (23) 


0{zi  -  g{zi))-^(j)o,o,i{zi)  +  A(1  -  g{zi))  +  C  (^1  -  ^  g(^i)) )  '^0’0,i(2i)  =  0.  (24) 


In  order  to  solve  the  differential  equation  (24),  we  first  examine  the  function  k{zi)  =  zi  —  g{zi)  = 
zi  —  Q{zi,  g{zi)).  Note  that 


d  _  d  _  a^(l-6*(0)  _  fe*(U  ^ 

c^2:i  dzi  aC- A(1 -6*(0)(a  +  0  5%  “ /» 

where  p  is  defined  in  Eq.  (1).  We  then  observe  that,  for  />  <  1,  the  quantity  k{zi)  never  becomes 

zero  in  \zi\  <  1,  while  for  p  >  1  this  quantity  has  one  and  only  one  zero,  call  it  /3,  such  that 

f3  G  (0, 1)  (see  for  example  [12]  or  [18]).  Now  rearranging  Eq.  (24),  define  the  function 


h{zi)  = 


A(1  -  5(^1))  +  e  (1  -  ^+A(l-g(2i))) 


2i  -  g{zi) 


and  note  that,  for  p  <  1, 

1-  u  AC(l-y(0)(a  +  0  _  Cp  ^ 

Thus,  we  conclude  that  h{zi)  is  analytic  on  the  open  disk  jzij  <  1,  and  h{zi)  can  be  defined  at  the 
point  zi  =  1.  Therefore,  on  the  closed  disk  jzij  <  1,  the  differential  equation. 


^<('0,0,1(21)  +  h{zi)(j)ofi,iizi)  =  0, 


is  verified  by  the  function 


<('0,0,1(21)  =  iLexp  -- 


1  ri  A(1  -  g{u))  +  C  (1  -  ^+A(r-gW)) 


g{u)  -  u 
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where  K  is  a  constant  of  integration.  For  p  <  1,  1(2^1)  in  Eq.  (25)  makes  all  generating  functions 

analytic  in  \zi\  <  l,\z2\  <  1.  In  particular,  (2),  (3),  and  (4)  are  obtained  up  to  the  multiplicative 
constant  K.  To  obtain  this  constant,  we  apply  the  normalization  condition,  +  V’o,o(l)  1)  + 

■01^(1, 1)  =  1,  which  yields 


«gb*(0-A(l-fe*(0)(a  +  0 


'/’o, 0,1(1), 


the  steady  state  probability  that  the  server  is  not  failed  and  idle. 

To  see  that  p  <  1  is  also  necessary  for  system  stability,  assume  that  p  >  1  and  the  system  is 
stable.  By  Lemma  1,  when  p  >  I,  the  function  zi  —  Q{zi,  g{zi))  has  one  and  only  one  zero,  denoted 
by  /?,  with  0  <  /?  <  1.  Thus,  the  function  h{zi)  is  not  analytic  in  l^ij  <  1.  Substituting  zi  =  P  in 
Eq.  (24)  yields 


A(l-/3)+dl- 


a 


0o,o,i(/3)  —  0, 


a  +  A(1  —  P) 

which  implies  that  0o,o,i(/3)  =  0  since  the  coefficient  of  0o,o,i(/3)  does  not  equal  zero.  We  must 
therefore  conclude  that 


0O,O,l(/3)  —  P^T^OfiJ,!  —  0. 

j=0 


However,  it  is  not  possible  to  find  positive  values  7ro,o,j,i,  j  >  0,  to  satisfy  the  above  relation,  so  the 
system  is  not  stable;  thus,  a  contradiction. 

Einally,  suppose  that  p  =  1  and  the  system  is  stable.  When  p  =  1,  the  function  zi  —  Q{zi,  g{zi)) 
never  becomes  zero  in  the  unit  disk  \zi\  <  1;  however,  in  this  case 


lim  h(zi)  =  — - - =  —00. 

2:1^1  1  ~  P 


Differentiating  Eq.  (24)  and  evaluating  at  the  point  zi  =  1,  we  obtain  the  relation 


00, 0,1(1)  —  0, 

21=1 


implying  that  0o,o,i(l)  =  0  since  ^g{zi)\^^^-^  =  1  when  p  =  1.  However,  this  contradicts  the 
hypothesis  that  the  system  is  stable.  Therefore,  we  conclude  that  the  system  cannot  be  stable 
unless  p  <  1.  m 
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Using  Theorem  1,  we  next  obtain  the  joint  distribution  of  the  orbit  and  normal  queue  size,  as 
well  as  the  distribution  of  the  overall  system  size,  independent  of  server  status. 

Corollary  1  For  p  <1,  the  probability  generating  funetions  G{zi,Z2)  and  H{z)  are  given  by 


^  J  X{1- B{z2))[a  + ^zi+ X{1- Z2)]{Z2- zi){l- g{zi))[a  +  ^  + X{1- g{zi))] 

L.(Zi,Z2)  —  <  - p - ; - : - ^ 

{  [aizi)  -  2;i][a  + A(1  -fl(2i))]  I (22  -  B[z2))[a  +  A(1  -  22)]^  + A(1  -  22)]  -  a^2i(l  -  B(22))| 

A(1  —  B{Z2y)\oi  +  ^2l  +  A(1  —  22)](1  —  22)[q;  +  C  +  -^(1  “  Z2)]  ^  [q  +  C  +  A(1  —  22)]  ] 

[a  +  A(1  -  22)]  |(22  -  B{z2))\a  +  A(1  -  Z2)][^  +  A(1  -  22)]  -  a?2i(l  -  B(22))|  “  ^2)]  j 


X  <^0,0.1(21)  (26) 


UjIIjUj 

[a  +  ^  +  A(1  -  z)]  [a^B{z)  +  \B{z){l  -  z)[a  +  C  +  A(1  -  z)]) 

H{z)  =  - ^ - ^4,o,i(^).  (27) 

[a  +  A(1  —  z)]  la^B{z)  —  X{z  —  B{z))[a  +  .^  +  A(1  —  z)]  | 

Using  standard  methods,  Eqs.  (26)  and  (27)  can  be  used  to  obtain  the  mth  moment  {m  >  1)  of  R, 
Q,  and  N,  respectively,  as  well  as  their  probability  distributions.  The  first  moments  are  provided 
in  the  following  corollary. 

Corollary  2  The  steady  state  mean  orbit  size,  mean  normal  queue  size,  and  mean  number  in 
system  are  respeetively  given  by 

ey(0K-A(l-r(0)]  +  («  +  0[A(l-y(g))-g^1  d  .28) 
^  ^  l-p  a  +  C  6*(0[«e-A(l-6*(0)(a  +  0]  0  ^  ^ 


E(Q)  =  A 


C^b*{C)  -  (g  +  C?[CB'  -  A(1  -  b*(0)] 
C6*(0(«  +  0  [a^  -  A(1  -  6*(0)(«  +  0]  ’ 


^  ^  C6*(0(«  +  0[«C&*(e)-A(i-6*(0)(«  +  0]  e{i-py  ^  > 


where 


d  ~  r°° 

B'  =  - — B{z2)  =  a  /  xe~^^b{x)da 

dz2  ^,=1  Jo 


Let  S  denote  the  duration  of  an  arbitrary  service  time  with  c.d.f.  B.  Then  as  ^  ^  0  in  Eqs.  (28) 
to  (30),  it  can  be  shown  that 


E{R)  =  0, 
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and 


A2e(52) 

2(1  -  AE(5))’ 


E(g)  = 


+  W^y 

These  expressions  are  consistent  with  results  for  the  standard  M/G/1  queue  with  no  failures  and 
no  retrials.  Furthermore,  denote  by  Wr,  Wq  and  W,  the  time  spent  in  the  orbit,  normal  queue 
and  system  by  an  arbitrary  customer  in  the  long-run,  respectively.  The  expected  values  of  these 
random  variables  are  obtained  by  applying  Little’s  Law  to  Eqs.  (28)  to  (30),  i.e.,  E(VFr)  =  A“^E(i?), 
E{Wq)  =  A-^E(Q),  and  E(1T)  =  X-^E{N). 

Corollary  2  also  confirms  that  p  <  1  is  necessary  for  the  stability  of  R  (and  N),  and  by  (29) 
we  see  that  pi  <  1  is  necessary  for  the  stability  of  Q.  That  is,  the  normal  queue  can  be  stable 
even  if  the  orbit  stability  condition  is  violated.  Owing  to  the  nature  of  the  orbit  dynamics,  retrial 
customers  are  subordinate  to  normal  customers  and  may  be  served  only  when  the  server  is  idle 
and  operational.  Hence,  normal  queue  customers  experience  a  greater  effective  service  rate  than 
do  retrial  customers,  and  thus,  it  is  possible  that  the  orbit  may  continue  to  grow  while  the  normal 
queue  remains  stable. 

Finally,  we  characterize  the  steady  state  distribution  of  the  server’s  status  by  directly  applying 
the  results  of  Theorem  1.  Let  p/,  pp,  and  pB  respectively  denote  the  limiting  probability  that  the 
server  is  idle,  failed,  or  busy. 


Corollary  3  For  p  <  1,  the  steady  state  distribution  of  the  server’s  status  is  given  by 

X{l-b*{0) 


a 


PI  =  lim  = 

zi^i  a  +  f 


PF  =  hm  V’o,o(^i,^2)  =  , 

21^1  a  +  c 

22^1 


and 


PB  =  lim  Tpipizi,  Z2)  = 

Zl^l 
22^1 

In  the  following  section  we  show  that  the  orbit  and  system  size  can  be  stochastically  decomposed 


before  considering  the  optimal  selection  of  the  retrial  and  repair  rates  in  section  5. 


A(l-y(0) 

Cb*{C) 
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4  Stochastic  Decomposition 


In  this  section,  we  demonstrate  that  both  the  orbit  and  system  size  exhibit  a  stochastic  de¬ 
composition  property  which  has  been  observed  for  the  system  size  distribution  of  many  M/G/1 
models  including  those  with  vacations,  retrial  queues,  and  breakdowns  (cf.  [4,  11,  16,  27]).  Falin 
and  Templeton  [12]  provide  several  stochastic  decomposition  results,  including  the  decomposability 
of  the  vector  of  server  status  and  the  orbit  size,  in  the  standard  M/G/1  retrial  queue  (i.e.,  one  with 
no  infinite  waiting  space  and  no  server  breakdowns). 

Allowing  0  ^  oo  in  our  model  yields  a  model  in  which  retrial  customers  instantaneously  attempt 
to  re-access  the  server  (i.e.,  an  instantaneous  feedback  model).  Let  R  denote  the  steady  state  orbit 
size  in  the  instantaneous  feedback  model,  and  denote  the  generating  function  of  R  by  K{z^)  for 
\z\  <  1.  Let  N  denote  the  steady  state  total  number  of  customers  in  the  system  in  the  instantaneous 
feedback  model  and  denote  its  generating  function  by  K{z^)  for  \z\  <  1.  Finally,  let  F  be  a  random 


variable  whose  generating  function  is  given  by 


E{z  )  =  exp  <  -- 


1  cl  A(1  -g(ir))  +e(l  -  ^+;,(f_g(^))) 


g{u)  -  u 


du  >  ,  1^1  <  1. 


The  following  two  propositions  describe  the  decomposability  of  the  orbit  and  system  size  distribu¬ 


tions. 
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Now  since  =  G{zi,  1),  we  may  write 


^G(^;i)exp  < 

E{z{^)E{zY) 

E(zf+'^). 


1  ri  A(1  -g(n))  +^(1  -  (,))) 

si,  - WT-n - 


Similar  behavior  may  be  observed  for  the  steady  state  system  size  as  noted  in  Proposition  2. 

Proposition  2  The  random  variable  N  may  be  expressed  as  the  sum  of  two  independent  random 
variables,  one  of  whieh  is  the  steady  state  system  size  in  the  instantaneous  feedbaek  model  and  the 
other  is  V,  i.e., 

N  =  N  +  V.  (33) 

Proof.  The  proof  is  analogous  to  that  of  Proposition  1.  Note  that  by  setting  zi  =  Z2  =  z  in 
Eqs.  (2),  (3),  and  (4)  we  obtain 

E{z^)  =  H{z)  =  +  i>o,o{z,  z)  +  ztfi^i^z,  z). 

The  generating  function  of  N  is  given  by 

E(z^)  =  lim  H{z)  =  Ah{z) 

9^00 


Ah{z)  = 


where 

a(l  -  p)[a  +  i  +  A(1  -  z)\  ^afB{z)  +  \B{z){l  -  2:)[a  +  ^  +  A(1  -  2;)]| 
(a  +  0[«  +  -^(1  “  ^)]  -  A(2:  -  B{z))[a  +  ^  +  A(1  -  2:)]| 

The  generating  function  for  N  is  the  product  of  these  two;  that  is, 

E(.-,  = 

y  0  Jzi  9{u)  -  u  j 

=  E(z^)E(z'") 

=  E(z^+'^). 
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In  the  next  section  we  formulate  an  optimization  problem  for  the  selection  of  the  optimal 
retrial  and  repair  rates  that  minimize  the  long-run  average  cost  of  operating  the  queueing  system, 
subject  to  a  budget  constraint.  We  also  provide  two  illustrative  examples  using  distinct  service 
time  distributions. 

5  Optimal  Retrial  and  Repair  Rates 

We  now  consider  the  simultaneous  optimal  selection  of  the  retrial  and  repair  rates  that  minimize 
the  long-run  average  operating  costs.  The  cost  function  includes  the  cost  of  service,  the  cost  of 
holding  customers  in  the  normal  queue,  and  the  cost  of  holding  customers  in  the  orbit.  Required 
for  the  optimization  are  the  queueing  performance  measures  IE(i?),  IE((5),  E(lTi{),  and  E(1Tq),  as 
well  as  the  expected  number  of  customers  in  service,  A(l  —  6*(^))/^6*(^),  and  the  expected  time  to 
complete  service,  (I  —  b* {^)) / ^b* {^) .  The  cost  per  unit  time  per  customer  in  service  is  cs  while  the 
holding  costs  per  unit  time  per  customer  in  the  orbit  and  normal  queue  are  respectively  denoted 
by  cr  and  cq.  The  coefficient  cg  is  the  cost  of  one  “unit”  of  retrial  rate  while  Cq  is  the  cost  of  one 
“unit”  of  repair  rate.  Using  Eqs.  (28)  and  (29),  we  solve  the  optimization  problem 


Minimize 
Subject  to 


«)  =  +  crE{R)E{Wr)  +  cqE{Q)E{Wq) 

A(l-6*(0)(a  +  0-«?&*(e)  <  0 

CgO  -|-  CaCX  <  D 
0,a  >0 


(34) 

(35) 

(36) 


where  D  is  a  fixed  budget  {D  <  oo).  Constraint  (34)  enforces  the  stability  condition  (p  <  1) 
and  (35)  is  a  budget  constraint  that  limits  the  attainable  repair  capacity  and  the  rate  at  which 
interrupted  customers  may  attempt  to  re-access  the  server. 

5.1  Convexity  Analysis 

The  uniqueness  of  a  global  solution  to  the  above  optimization  problem  can  be  established  by 
showing  it  is  a  convex  program.  The  existence  of  a  solution  will  then  be  shown  directly. 
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The  feasible  region,  defined  by 


X  =  {(0,  a)  :  A(1  —  b*{^)){a  +  C)  “  <  0;  cgO  +  Cq-ck  <  D;  9,a  >  0}, 

is  a  convex  set  since  it  is  defined  by  a  finite  set  of  linear  constraints.  Strict  convexity  of  the  objective 
function,  which  will  now  be  proved,  will  complete  the  uniqueness  proof.  Expanding  the  terms  of 
C{6,  a)  gives 


eb*{C)  -  {a  +  C)^[CB'  -  A(1  -  6*(0)] 


^25* (^)2 

+  c^A  faA(l-6*(0) 


+  0  [«?  -  A(1  -  b*{0){a  +  0]  J 

Cb*m  -  A(1  -  b*m  +  (g  +  g)[A(l  -  b*{0)  - 
b*{0H  -  A(1  -  b*{()){a  +  -  A(1  -  b*{0){a  +  0] 


+ 


AC(«  +  g)(l-fc*(0) 

-A(l-6*(0)(a  +  0] 


(37) 


The  first  term  on  the  right-hand  side  (r.h.s.)  of  Eq.  (37)  depends  on  neither  9  nor  a,  and  hence, 
does  not  affect  the  convexity  of  C.  To  prove  the  convexity  of  the  other  two  terms,  the  following 
functions  are  defined: 


f  f .  ^b*{m  -  A(i  -  b*m  +  (« + o[A(i  -  y (g))  - 

>  [ae  -  A(1  -  b*{0){a  +  0][«e6*(e)  -  A(1  -  6*(0)(«  +  0]  ’ 


/-R2(^)«) 

fnio,  a) 


Ag(«  +  0(l-b*(0) 

0K&*(O-a(i-6*(O)(«  +  O]’ 


A(l-y(0) 

6*(0 


/iJi(a)  + 


fb*{^)-{a  +  Cf[C&-X{l-b*{m 

(a  -h  0  [a^  -  A(1  -  6*(0)(a  +  0] 


The  next  two  lemmas  are  needed  to  prove  the  strict  convexity  of  C{9,  a)  on  X. 


Lemma  2  The  function  f‘^{9,a)  is  strictly  convex  on  X. 


Proof.  We  first  establish  the  positivity  of  the  quantity 


X{l-b*{C))-CB' 


A  1- 


{f,x  +  l)b{x)e 
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Since  its  derivative  with  respect  to 


/  x^b{x)e  ^^dx, 

Jo 

is  strictly  positive,  A(1  —  b*{^))  —  is  strictly  increasing  for  G  [0,oo)  and  thus  attains  its 
minimum  value  of  zero  at  the  left  endpoint  .^  =  0.  Now  the  second  derivative  of  with  respect 
to  a  is  given  by 


_  y  f  -  A(1  -  -  A(1  -  ^(0)] 

da^  ^  \  [a^b*(^)  -  A(1  -  b*(^))(a  +  C)]^[a^  -  A(1  -  b*(0)(a  +  0]^ 

[A(l  -  b*(0)  -  ^B']  {«[^  -  A(1  -  b*{m[^b*{0  -  A(1  -  y (0)]  +  A2g(l  -  y (g))^} 
[a^b*{0  -  A(1  -  b*{C)){a  +  0]^K  “  -^(1  “  &*(?))(«  +  0]^ 

A(1  -  b*{0){a  +  -  A(1  -  6--(0)][A(l  -  b*{0)  -  mHb*{0  -  A(1  -  b*{0){a  + 

Hb*{0  -  A(1  -  6*(0)(a  +  -  A(1  -  6*(e))(«  +  0]^ 

A^r (0(1  -  y (g))[g  -  A(1  -  b*{0)AHb*{0  -  A(1  -  b*{0){a  +  Q]^ 

[«e&*(0  -  A(1  -  6*(0)(a  +  -  A(1  -  6*(0)(«  +  0]^ 

A(1  -  b*ma  +  mb*{0  -  A(1  -  fen0)][A(l  -  y (0)  -  -  A(1  -  y (0)(«  +  Q]^ 

Hb*{0  -  A(1  -  6*(0)(a  +  -  A(1  -  6*(0)(«  +  0]^ 

Acr(0(i  -  b*m(  -  A(i  -  b*{omb*{o  -  A(i  -  b*mH  -  a(i  -  b*mc^ + ) 
[a^b*{C)  -  A(1  -  6*(0)(a  +  0]3[ae  -  A(1  -  b*{0){a  +  0]3  /  ' 

It  follows  from  the  stability  condition  that  [a^  — A(1  —  6*(0)(a  +  0]  >  0,  [C&*(0  —  A(1  —  6*(0)]  >  0, 

and  —  A(1  —  6*(0)]  >  0.  Since  A(1  —  6*(0)  —  ^B'  is  nonnegative, 


and  thus,  fn^ia)  is  strictly  convex  for  all  a  >  0.  Taking  second  partial  derivatives  of  with 

respect  to  9  and  a,  respectively,  and  assuming  stability,  shows  that 

d^fR,{e,a)  _  2A0i-&00)(«  +  0 

dO"^  03[a^6*(O  —  A(1  —  6*(0)(a  +  0] 


and 

d^fRAe,a)  _  2Xeb*A){l  -  b*m^b*{0  -  A(1  -  b%m 
da^  0[a^6*(O  —  A(1  —  6*(0)(a  +  0]^ 


Let  H(0,  a)  denote  the  Hessian  matrix  of  fR2{9,  a).  Then  it  can  be  shown  that 


det(H(0,  a)) 


A^0b*(0(l  -  b*{0A  {30b*(0  +  4[aCb*{C)  -  A(1  -  b*{0){a  +  Q]} 
(94[a^6*(0  -  A(1  -  6*(0)(a  +  p]^ 


Hence,  fR^{6,a)  is  strictly  convex  on  X.  Consequently,  fR  =  /r^  +  /r^  is  strictly  convex  on  X, 
and  thus,  is  strictly  convex  on  X.  ■ 


18 


Lemma  3  The  function  fqia)  is  strictly  convex  for  all  a  >  0. 


Proof.  The  second  derivative  of  fq  with  respect  to  a  is  given  by 

^  2^2  /  (« + OK  -  A(i  -  b*m{eb*{o  -  (« + mcB’  -  a(i  -  b*m} 
da^  ^  \  {a  +  Cna^-X{l-b*{0){a  +  CW 

-  A(1  -  6*(0)(«  +  e)]{2e'  +  3[a^  -  A(1  -  b*{^)){a  +  0]}  1 

(«  +  03K-A(l-6*(0)(a  +  CF  /■ 

It  follows  from  the  stability  condition  that  [a^  —  A(1  —  +  >  0  —  A(1  —  6*(0)]  >  0, 

and  Lemma  2  ensures  that  A(1  —  b*{f))  —  ^B'  is  nonnegative.  Hence, 

da^  ’ 

which  establishes  strict  convexity  of  fq  for  a  >  0.  The  strict  convexity  of  /^(a)  follows  directly.  ■ 
The  following  theorem  is  the  main  result  of  this  section. 

Theorem  2  The  cost  function  C{6,a)  is  strictly  convex  on  X. 


Proof.  The  proof  follows  directly  from  Lemmas  2  and  3.  In  particular,  the  strict  convexity  of 
/g(a)  for  all  a  >  0  ensures  that  f^{0,  a)  +  fq{a),  and  thus  C{6,  a),  is  strictly  convex  on  X.  m 
Theorem  2,  along  with  the  convexity  of  X,  show  that  the  optimization  problem  is  a  convex 
program  (CP).  A  CP  guarantees  that  any  stationary  point  (Karush-Kuhn- Tucker  point)  is  a  global 
minimizer,  but  to  ensure  existence  of  a  solution,  the  feasible  region  must  be  closed  and  bounded. 
Although  X  is  bounded,  it  is  not  closed.  To  circumvent  this  complication,  we  note  that,  for  all 
values  of  6  and  a, 


dC{9,  a)  _  -2cRW{a  +  ^(1  -  b*{i)) 

de  e^[aib*{i)-\{l-h*{i)){a  +  e)] 
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and 


dC{e,a)  _ 

da 

-2cR\^i{i-h*{i))fR{e,a) 


+ 


eb*{i) 


A(1  -  b*{0)  |(«  +  6[A(1  -  h*m  -  ^B']  +  ih*{m  -  A(1  -  6=* 
h*{i)[aih*{i)  -  A(1  -  6*(0)(a  +  -  A(1  -  6*(0)(a  +  0] 

-  6*(0)  -  -  A(1  -  6*(0)(«  +  0]"} 


a 


+ 


e[aib*{i)  -  A(1  -  6*(0)(«  +  e)]  b*{i)[aib*{i)  -  A(1  -  b*{i)){a  +  0][aC  -  A(1  -  6*(0)(a  +  0]" 

' eb*{0  -  (a  +  -  A(1  -  b*m  +  2e^*(0[«e  -  A(1  -  6*(0)(«  +  6] ' 


-  2AcQ/g(a) 


6*(C)2(a  +  -  A(1  -  b*(0)(a  +  0]^ 


<  0. 


Hence,  the  cost  function  C(0,  a)  is  monotonically  decreasing  in  both  6  and  a  and  is  bounded  below 
by 


A(1  -  b*{C)f  .  A3(l  -  6*(0)'[A(1  -  b*{0)  -  X[CB'  -  (1  -  b*{m^ 

_  6*(^))]2[^  _  A(1  -  6*(e))]2+"Q^25*(^)2[^  _  ^(1  _  b*m^’ 


Therefore,  the  budget  constraint  (35)  is  always  binding,  and  we  may  substitute 

Q _  B  CqQ; 

ce 


into  C{6,a).  Differentiating  the  objective  function  with  respect  to  a  and  setting  it  equal  to  zero, 
the  optimal  repair  rate,  denoted  by  a*,  is  the  unique  root  that  satisfies  Eq.  (34).  Subsequently, 
the  equation, 

^  D  -  Cga* 

Ce 

may  be  solved  to  obtain  the  optimal  retrial  rate  B* .  Thus,  (0*,  a*)  is  a  stationary  point  of  C  on  X, 
and  hence,  the  global  minimizer. 


5.2  Numerical  Examples 


We  now  illustrate  the  solution  procedure  in  two  distinct  scenarios.  In  the  first  case,  we  assume 
the  service  times  are  exponentially  distributed  with  positive  rate  parameter  /x.  In  this  case,  the 
Laplace  transform  of  the  service  time  distribution  is  given  by 
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which  gives 


and 


r(0  = 


B'  = 


(/^  +  0 

Xfj, 


In  the  second  scenario,  we  assume  the  service  times  are  uniformly  distributed  on  the  interval 
(0,2//i).  In  this  case,  the  Laplace  transform  of  the  service  time  distribution  is 

/r(l  - 


b*{s)  = 


2s 


so  that 


and 


b*{0  = 


B'  = 


Xfi 


1  —  exp 


/i(l  - 

2e 

-2C 


2£ 

1  +  — 


The  cost  coefficients  in  both  cases  are:  cg  =  Ca  =  cs  =  I,  and  cr  =  cq  =  100.  The  remaining 
parameters,  as  well  as  the  optimal  solutions  (indicated  by  *),  are  specified  in  Table  1. 


Table  1:  Optimal  repair  and  retrial  rates  for  two  numerical  examples. 


Service  time  distribution 

A 

D 

9* 

a* 

C{9*,a*) 

Exponential  (/x) 

5.0 

10.0 

1.5 

30.0 

11.23 

18.77 

5.85 

Uniform  on  (0,2//i) 

5.0 

10.0 

1.5 

30.0 

11.69 

18.31 

5.18 

As  expected.  Table  1  verifies  that  the  budget  constraint  is  binding  at  the  optimal  solution.  In 
particular,  with  ce  =  =  1,  we  have  9*  +  a*  =  D  =  30.0  in  both  examples. 
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